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PART  A 


Molecular  Dynamics  and  Ionic  Associations  of  LiAs  F6 
in  4-Bulyrolactone  Mixtures  with  2-Methyl  Tetrahydrofuran 

Yoshifumi  Harada,1  Mark  Salomon 2  and 
Sergio  Petrucc »' 

Polytechnic  Institute  of  New  York 
Route  110 

Farmlngdale,  New  York  11735 

U.  S.  Army  ET  &  DL 
Power  Sources  Division 
Ft.  Monmouth,  New  Jersey  07703-5302 

ABSTRACT 

Solutions  of  LLAsF6  at  2fi8.2  K  were  studied  by  audiofrequency 
conductance  in  the  concentrations  range  10_4-10-2  mol  dm'3,  and 
by  radiofrequency  ultrasonic  absorption  in  the  concentration  range 
0.05  -  0.5  mol  dm-3.  The  solvents  employed  were  mixtures  of  4- 
butyrolactone  with  2-methyl  tetrahydrofuran  varying  In  composi¬ 
tions  from  mole  fractions  xbl  ~  0.10  to  Xbl  =  0.75.  In  dilute 
solutions  the  audiofrequency  conductivity  data  yield  ion  associa¬ 
tion  constants  which  appear  to  represent  both  contact  and 
solvent-separated  ion  pairs  for  mixtures  up  to  X^  =  0.36.  At 
higher  LLAsFs  concentrations  and  In  solvent  mixtures  of  composi¬ 
tion  XBL  =  0.75,  and  XgL  =  0.35,  the  ultrasonic  spectrum  shows 

1  Polytechnic  Institute  of  New  Yorlt,  Pvmlnfdale,  New  Yorlt. 
o 

U.  8.  Army  ET  ft  DL,  Ft.  Monmouth,  New  Jersey. 


3 


.o 


1 


a  single  relaxation  process  which  again  is  attributed  to  the  formar 
tlon  of  both  contact  and  solvent  separated  Ion  pairs.  For  = 
0.10,  the  ultrasonic  spectrum  Is  the  sum  of  two  Debye  relaxation 


processes 


identified 


equilibria 


2LIAsFfi  £+  LiAsF6...LiAsF6  ^  (L1AsF6)2. 


Introduction 


Over  the  past  five  years  L1AsF6  has  emerged  as  the  preferred  electrolyte 
for  use  in  rechargeable  lithium  batteries  [l].  A  major  reason  for  this  selection 
Is  that  L1AsF6  solutions  generally  show  much  higher  electrolytic  conductivities 
than  corresponding  solutions  of  other  lithium  salts  such  as  LICIO.,  or  LIBF4.  In 
addition  to  a  smaller  tendency  to  form  ion  pairs  (2-6),  LlAsF6  also  has  a 
smaller  tendency  to  form  dimers  (i.e.,  quadrupoles)  than  do  LICIO^  and  LiBF4 
(2,3).  The  use  of  ether  based  solvents  is  also  of  practical  interest  since  lithium 
anodes  appears  to  cycle  well  in  these  solvents  (l).  The  use  of  binary  mixtures 
of  an  ether  with  a  high  dielectric  constant  solvent  such  as  4-butyrolactone  (4- 
BL)  and  propylene  carbonate  (PC)  has  still  greater  practical  applications  since  it 
has  been  shown  (4,6,7  and  references  cited  therein]  that  electrolytic  conduc¬ 
tances  are  generally  greater  In  specific  mixtures  than  they  are  in  any  of  the  pure 
solvents.  The  present  study  Is  concerned  with  the  properties  of  LIAsF6  in 
binary  mixtures  of  4-BL  with  2-methyl  tetrahydrofuran  (2MeTHF).  The  sys¬ 
tems  were  studied  by  classical  audiofrequency  conductance  measurements  at 
low  concentrations  and  by  more  recent  relaxation  kinetic  methods  which  have 
the  advantage  of  being  able  to  isolate  (in  the  frequency  spectrum)  one  or  more 
multistep  processes  (8).  Parallel  studies  of  LLAsF6  in  pure  2MeTHF  and  pure 
4-BL  have  been  reported  earlier  (3,5]. 

Experimental 

Chemicals  LIAsF6  (USS  Agri-Chemicals  "Electrochemical  Grade")  was  dried  at 
60-  70*C  under  vacuum  for  24h.  2MeTHF  was  distilled  from  benzophenone  at 
reduced  pressure  (3)  or  was  distilled  from  type  4A  molecular  sieves  under  an 
atmosphere  of  pure  argon.  4-BL  (Aldrich  "Gold  Label"  product)  was  distilled 
from  4A  and  5A  molecular  sieves  at  reduced  pressure  as  described  In  (5,6].  All 
solvent  mixtures  were  prepared  by  weight. 


-  5  - 


Solvent  Properties  The  static  dielectric  constants  of  the  solvent  mixtures  were 
measured  at  25±0.1*C  by  the  comparison  method  [4-6],  and  a  precision  of 
0.5%  or  better  was  attained.  Densities  were  measured  at  25.0±0.1*C  with  a 
PAAR  DMA  45  digital  density  meter.  Viscosities  were  measured  at 
25.00±  0.02  *C  with  a  Lauda  precision  viscometer  with  a  precision  of  ±  0.1%. 
The  physical  properties  of  these  solvents  are  given  in  Table  1.  The  data  in 
Table  1  are  fitted  by  least  squares  to  the  following  smoothing  equations 

e(D)  =  6.137  +  26.222Xjjl  +  7.8952X|L  +  1.6939X|L,  (l) 

»?(cP)  =  0.4653  +  0.5691XBL  ~  0.04594X|l  +  0.7387X|L,  (2) 

d(gcm“3)  =  0.8475  +  0.2209XBL  +  0.06028X|l  -  0.005223X|L,  (3) 

Note  that  the  solvent  compositions  are  mass  %  in  Table  l,  and  the  smoothing 
equations  (1-3)  are  based  upon 'mole  fraction  of  4-BL. 

Audiofrequency  Measurements  Electrolyte  solutions  were  prepared  by  weight 
in  a  dry  box  and  conversions  to  volume  concentration  units  were  calculated 
from  measured  densities.  Aliquots  of  a  given  stock  solution  were  placed  in 
three  10  ml  flasks  which  were  then  used  for  the  conductivity  runs  using  succes¬ 
sive  weight  dilutions.  Three  Jones  type  conductivity  cells  were  used  in  combi¬ 
nation  with  an  Altex  RC-18  conductivity  bridge.  Electrolytic  conductance  were 
measured  at25.00*C  (accuracy  ±0.02K  )  is  described  in  [4-8]. 

Ultrasonic  Relaxation  The  equipment  and  procedures  have  been  described 
elsewhere  [3],  A  new  resonator  ultrasonic  cell  with  1  inch  diameter  quartz  cry¬ 
stals  was  used  as  it  allowed  measurements  over  a  broader  frequency  range 
(0.5-16  MHz),  and  a  better  signal  to  noise  ratio  than  the  cell  used  previously 
[3].  After  thermostating  it  at  25.00  *C  (precision  ±0.01K),  for  2  hours,  fre¬ 
quency  drifts  in  the  resonating  frequencies  stopped.  Solutions  were  kept  in 


desiccators,  used  within  24  hours  of  preparation,  and  were  exposed  to  the 
atmosphere  during  filling  the  resonator  cell  for  20-30  seconds  at  most. 

Results  and  Discussion 

The  molar  conductivities,  A(Scm2mol_1),  were  calculated  from  the  experi¬ 
mental  electrolytic  conductances  after  correcting  for  the  solvent  conductance, 
and  the  data  are  given  in  Table  2.  The  data  in  Table  2  were  analyzed  by  means 
of  the  Fuoss-Hsia  equation  (0)  using  the  method  proposed  by  Justice  [10]  and 
the  constants  derived  by  Fernandez-Prlnl  (llj, 

A  =  A00  -  S( arc) 1/2  +  Eac  flmc  +  J^Rjac  -  J2(R2)(ac)3/2  -  Aacy2  Kfl(4) 

where  all  terms  have  their  usual  significance.  The  thermodynamic  ion  associa¬ 
tion  constant  and  mean  molar  activity  coefficient  are  defined  In  the  usual 
manner: 

K£°=  (l-or)/ar2cy£  ,  (5) 

and 

flay*  =  -A(ac)1/2/(i  +  BRy(ac)1/2]  (6) 

The  distance  parameters  R,  and  Ry  in  Eqs.  (4)  and  (6)  were  set  equal  to  the 
BJerrum  distance  q  as  suggested  by  Justice  [10],  and  Eqs.  (4-6)  solved  by  a 
least  squares  method  similar  to  that  described  by  Kay  [12].  Values  of  q 
(expressed  in  A)  are  given  In  Table  3.  Table  3  also  gives  the  results  of  fitting 
the  parameters  A00  ,K£°,  and  R2  to  Eqs.  (4-6).  a  A  In  this  table  is  the  standard 
error  of  estimate  for  the  molar  conductivities.  Based  on  the  a  A  values,  the 
uncertainties  in  A  are  around  0.2% (96%  level  of  confidence). 

For  electrostatic  Interactions,  K“  can  be  calculated  from  the  BJerrum  equa- 


Kf°  - - —  /  R2exp[z2e2/ckTR]dR  , 

a  1000  \ 


(7) 


where  the  Integration  Is  carried  out  from  the  distance  of  closest  approach,  a,  to 
the  BJerrum  distance,  q.  For  a  contact  ion  pair,  the  distance  of  closest  approach 
can  be  simply  set  equal  to  the  sum  of  the  ionic  radii.  For  LiAsF6  a=  r+  +  r_ 
—4.44A,  and  in  which  case  the  calculated  I<£°  values  for  the  mixed  solvents  are 
all  much  larger  than  the  experimental  K^°  values  given  in  Table  3.  Calculated 

values  can  be  reduced  by  considering  solvent-separated  ion  pairs,  and  in 
which  case  a  =  r+  +  r_  +  d  where  d  is  the  diameter  of  the  solvent  molecule 
(4-7).  However,  taking  d^S^A  (the  approximate  solvent  diameter  of 
2MeTHF)  results  In  calculated  K£°  much  smaller  than  the  experimental  values. 
The  experimental  K£°  values  in  the  mixed  solvents  can  be  reproduced  from  Eq. 
(7)  by  taking  d~1.6A  In  fact,  K£°  values  for  LLAsF6  In  dim ethoxye thane  mix¬ 
tures  with  PC  or  4-BL  can  also  be  reproduced  within  experimental  error  by  tak¬ 
ing  d«l.GA(0).  This  low  value  of  1.6A  suggests  that  both  contact  and 
solve nt-separated  ion  pairs  exist  in  these  mixed  solvents  as  proposed  earlier  [6]. 
Additional  evidence  for  the  simultaneous  existence  of  both  types  of  ion  pairs 
comes  from  our  relaxation  studies  discussed  below. 

Turning  to  the  relaxation  studies.  Fig.  1  shows  a  representative  plot  of  a/f2 
as  a  function  of  frequency  f  for  L!AsF6  in  a  X^  =  0.75  mixture  at  25‘C.  The 
solid  line  is  the  calculated  Debye  functions  for  a  single  relaxation  process: 


a/f2  = 


l+(f/fr)5 


+  B 


(8) 


where  a  is  the  absorption  coefficient  of  sound  in  Npcm-1,  fr  the  relaxation  fre¬ 
quency,  and  A  and  B  are  fitting  parameters  (8,13).  The  data  and  parameters 
for  L1AsF6  In  the  XBL  =  0.75  mixture  are  given  In  Table  4.  The  relaxation 


parameter  A  In  Table  4  Is  quite  small  and  for  O.lmol  dm-3  LlAsF6  in  pure  4- 
BL  the  relaxation  effect  is  not  detectable.  This  is  consistent  with  the  conduc¬ 
tivity  results  (5)  which  show  that  LlAsFg  is  completely  ionized  in  4-BL  at  25*C. 
The  observed  single  relaxation  process  is  interpreted  in  terms  of  Jon  association 

kf 

Li+  +  AsFg  LLAsFg  ,  (0) 

for  which  one  can  write  the  rate  expression  (8) 

r" 1  =  kt6  +  kr  . 

In  Eq.  (10)  r-1  is  the  Inverse  of  the  relaxation  time  (r-1 
function  defined  by 

6  =  2'ycy.|  .  (ll) 

Values  of  0  were  estimated  from  Eqs.  (11)  and  (12)  taking 
K”  =  ao.SmorMin3  (obtained  by  interpolations  from  the  conductivity  data. 
Table  3).  In  the  Davies  equation  (14)  Ry  was  again  set  equal  to  the  BJerrum 
distance  q  (for  XBL  at  25*C,  q  =  9.08A). 


(10) 

=  27rfr),  and  0  is  a 


-AI*/2 

flny±  =  — .  -  +  Cl 


1+RyI1/2 


(12) 


Fitting  t_1  and  6  to  Eq.  (10)  by  the  method  of  least  squares  gives 
k,  =  (1.79  ±  0.49)  xl09mor1dm3s_1  and  kr  =  (2.85  ±  0.26)xl08s“ 1  with  a 
correlation  coefficient  r2  =  0.87.  These  rate  constants  give  an  ion  association 
constant  of  Kf°  =  k,/kr  =  6.3  ±  2.3  mol-1dm3  which  differs  from  that  of 
39.8mol_Idm3  obtained  from  the  conductivity  data.  Differences  of  this  magni¬ 
tude  are  typical  when  comparing  Kf°  values  from  conductivity  data  and  rate 


data,  and  are  probably  due  In  part  to  the  larger  experimental  errors  in  the  rate 
data  and  the  use  of  Eq.  (12)  for  calculating  activity  coefficients  in  highly  con¬ 
centrated  solutions. 


The  rate  constant  for  a  diffusion  controlled  Ionic  encounter  can  be  calcu 
lated  from  the  Debye-Smoluchowsky  equation  (15) 


,  8RT  b 

kD  = - * - —  ,  (13) 

D  30007/  i-e"b 

where  b  is  the  Bjerrum  parameter  |z+z_  |e2/qekT.  For  XBL  =  0.75  the  quanti¬ 
ties  c(30.85D)  and  t/(1.178cP)  were  obtained  from  Eqs.  (1-2),  and  thus 
kD  =  1.29X1010  mol-1dm3  s-1  which  Is  an  order  of  magnitude  greater  than  the 
experimental  value  of  1.70xl09morldm3s*1  for  k(.  This  result  leads  us  to  con¬ 
clude  that  for  L1AsF6  in  Xbl-0  .75  the  association  process  takes  place  with  a 
energy  barrier  larger  than  viscous  flow.  Elimination  of  solvent  from  the  first 
coordination  of  the  ions  is  the  likely  source  of  the  above  differences  in  the  k,'s. 
We  then  propose  the  Eigen  mechanism  (8]  for  the  association  process: 

ko  k, 

LI+  +  AsFg  +  S  Ll+*S'AsFg  j*  LlAsF6  +  S  ,  (14) 

k_0  k-l 

where  Li‘S*AsF6  is  the  solvent-separated  (intermediate)  ion  pair. 

Considering  LlAsF6  in  mixtures  where  XBL  =  0.35,  Fig.  2  shows  a 
representative  plot  of  the  excess  sound  absorption  per  wavelength 
H  =  (a-Bf2)u/f  =  aexcX  against  the  frequency  f.  Here  X  =  u/f  Is  the 
wavelength  and  B  is  the  value  of  o/f2  for  f»fr  (cf.  Eq.  (8)  ).  The  solid  line 
in  Fig.  2  Is  the  Debye  function  for  a  single  relaxation  process 


M  ~  2  M max 


t/tr 


(15) 


where  /imax  =  Aufr/2.  Table  5  gives  and  the  sound  velocity  u  for 


L1AsF6  solutions  in  XBL  =  0.35  mixtures  at  25*C.  From  Table  5  It  Is  seen  that 
within  experimental  error,  fr  is  independent  of  the  L1AsF6  concentration,  and 
that  ftmax  Is  a  linear  function  of  concentration.  From  the  conductivity  data,  K£° 
Interpolated  from  a  plot  of  log  K“  vs.  l/e  Is  457mor1dm3,  which  means  that 
the  concentration  of  free  ions  In  the  solutions  of  composition  Xql—  0.35  is 
small.  Under  these  conditions  the  scheme  presented  in  (14)  above  probably 
reduces  to 


LI*S-AsFa  ^  L1AsF6  +  S  ,  (16) 

k-t 

l.e.,  to  a  pseudo  first  order  process,  the  solvent  S  being  In  large  excess.  For 
this  pseudo  first  order  process,  the  rate  equation  Is  therefore  (8) 

r~ 1  =  kj  4-  k_!  =  (4.3  ±  0.2)xl0®s_1  . 

The  value  of  /imax  Is  also  given  by  (2) 


IT  (AV,)8  K, 

2/3,  RT  (1+K,)2  C  * 


(17) 


where  K£°  =  K^l+Kj)  =  457mol-1dm3  Kq  and  Kj  are  the  equilibrium  con¬ 
stants  for  the  two  steps  given  In  Eq.  (14)  and  the  adiabatic  compressibility 
/3S  =  (d0u2)-1  =  62.9x10“  12g"  lcm  s2.  At  the  present  time  K1  is  not  known 
and  AVa  cannot  therefore  be  calculated. 

Considering  L1AsF8  solutions  at  25*C  In  Xg^  =  0.10,  a  representative  plot 
of  y,  as  a  function  of  f  Is  shown  In  Fig.  3.  The  data  in  this  figure  cannot  be 
fitted  satisfactorily  in  terms  of  a  single  Debye  relaxation.  The  solid  line  In  Fig. 
3  was  obtained  by  fitting  the  data  to  two  Debye  functions  given  by  the  relation 
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^max.I 


f/fl 

l-Kf/fi)2 


l+(f/fll)2 


(18) 


Table  6  gives  the  parameters  fimaxi  »f|  *  A*  maxii  »fn»  the  background  absorption 
B,  and  the  velocity  u  as  a  function  of  LlAsFg  concentration.  The  inset  In  Fig.  3 
shows  the  tail  of  the  a/f2  vs.  f  corresponding  plot,  with  the  solid  line  calculated 
from  the  parameters  of  Table  6,  indicating  that  the  choice  of  B  is  rather  une¬ 
quivocal  and  subject  to  an  error  not  larger  than  ±0.5  to  ±  1.0x10"  17cm-1s2: 
i.e.,  the  second  relaxation  Is  required  to  accurately  describe  the  processes  In 
Xbl  =  0  .10  solutions.  From  the  plot  of  log  vs.  l/c  based  on  the  data  In 
Table  3,  for  XBL  =  0.10  we  find  that  Kf°x  =  2xl04mor  ‘dm3  which  means  that 
ion  association  is  virtually  complete,  and  that  no  ultrasonic  relaxation  can  arise 
from  dissociation  to  free  Ions.  Outer-sphere  Inner-sphere  processes  are  first 
order  (or  pseudo  first  order  as  described  above  for  LiAsF6  in  XBl  -  0.35  mix¬ 
tures)  which,  for  LiAsF6  in  XBL  =  0.10  mixtures,  cannot  explain  the  concentra¬ 
tion  dependencies  of  the  relaxation  frequencies  fj  and  f„.  In  order  to  explain 
these  concentration  dependencies  we  propose  (see  also  (3))  that  the  two  Debye 
processes  are  associated  with  the  dimerization  of  Ion  pairs  to  form  quadrupoles 
according  to: 


k3 

2LIAsF6  LLAsF6  •  •  •  S  •  •  •  LlAsFg  (LiAsF6)2  ,  (19) 

k_ 2  k_3 

where  the  dimerization  (or  quadrupole  formation)  constant  is  given  by: 


K,  =  K*(l+K3) 


(20) 


In  (19)  we  propose  a  solvent  separated  intermediate  quadrupole, 
LIAsF6...S...L1AsF6,  on  the  basis  of  the  two  observed  Debye  relaxations. 
LlAsFg  solutions  in  pure  2MeTHF  (3)  show  only  a  single  Debye  process,  and 


at  this  time  we  cannot  determine  which  solvent  component  S  represents  in.Eq. 
(10).  According  to  the  scheme  In  Eq.  (10),  one  would  expect  the  two  relaxa¬ 
tion  times  to  be  described  by  the  relation  (3) 

rui=  (S±  (Ss-4P)1/2)/2  (21) 

where 

S  =  rf 1  +  rfi 1  =  4k2[AB]  -f  k_2  +  k3  +  k_3  (22) 

and 

p  =  rf1  r,7 1  =  4k2[AB](k3  +  k_3)  +  k_2k_3  (23) 

When  the  dimerization  constant  Kq  is  small,  then  the  ion  pair  concentration 
[AB]  can  be  safely  approximated  by  the  stoichiometric  concentration 
c/mol  dm-3.  Linear  regression  of  S  vs.  c  gives  a  determination  coefficient 
r2  =  0.085,  a  slope  =  1.20x10®,  and  an  intercept  =  8.15x10®.  Linear  regres¬ 
sion  of  P  vs.  c  gives  r2  =  0.072,  slope  =  4.88xl017,  and  Intercept  = 

l.lOxlO17.  Similar  to  the  treatments  for  LiBF4  In  pure  2MeTHF  (3),  the 
above  slopes  and  intercepts  are  used  to  calculate  the  following  rate  constants 
and  equilibrium  constants: 

kj  =  3.2xl0®mor1dm3s~l;  k_2  =  4.4xl0®s" 1 
=  k2/k_2  =  0.73  mol-,dm3 
k3  =  l.lxl0®s_1;  k_3  =  2.7xl0®s_1 

K3  =  0.41 

K„  =  ^(l+Kg)  =  l.o  mol"  ‘dm3 

This  value  of  I<q  in  the  =  0.10  mixture  is  consistent  with  the  value  of 
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1.8  mol" ‘dm3  for  L1AsF6  in  pure  2MeTHF  (3). 

Conclusions 

The  combination  of  two  completely  different  methods  of  Investigation, 
albeit  In  different  concentration  ranges,  has  shown  a  consistent  description  of 
the  status  of  the  electrolyte,  as  being  practically  all  Ionized  (up  to  c«0.lM)  in 
4-BL,  and  practically  all  associated  with  a  small  percentage  of  dimerization,  in 
2MeTHF.  The  addition  of  4  BL  to  2MeTHF  causes  solvation  of  the  dimers  and 
Increases  the  ionization  of  the  electrolyte.  Perhaps,  the  potentially  chelating 
ability  of  4  BL  causes  preferential  solvation  of  Li+,  a  point  that  would  be  desir¬ 
able  to  check  experimentally  by  tools  others  than  the  ones  used  In  the  present 
study. 

It  Is  clear  from  the  above  Investigation,  that  no  method  has  the  capability 
of  giving  a  complete  picture  of  a  given  system.  Rather,  the  Information  gath¬ 
ered  by  several  parallel  investigations  by  different  technologies  will  give  enough 
complementary  data  as  to  offer  an  adequate  insight  for  a  molecular  description 
of  the  system. 
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Fig.  1  Representative  plot  of  the  quantity  (a/f2)  vs.  frequency  f,  for 


L1AsF6  in  the  solvent  mixture  2MeTHF-4-BL  of  composition 


X4R,  =  0.75  at  t  =  25*C. 


Fig.  2  Representative  plot  of  the  quantity  ft  vs.  f  for  LlAsF6  in  the  sol¬ 


vent  mixture  2MeTHF-4BL  of  composition  X4BL  =  0.35  at 


t  =  25*C. 


Representative  plot  of  the  quantity  /i  vs.  f  for  LiAsF6  in  the  sol¬ 


vent  mixture  2MeTHF-4BL  of  composition  X4BL  =  0.10  at  25*C. 
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Molecular  Relaxation  Dynamics  and  Ionic  Association 
of  LiBF4  in  1, 2-Dime thoxymethane 

by 
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Depts.  of  Chemistry  and  Electrical  Engineering 
Polytechnic  Institute  of  New  York 
Long  Island  Center 
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Abstract 

Electrical  conductance  data  for  L1BF4  in  1,2-Dimethoxymethane  (DMM) 
at  25°C  reveal  the  electrolyte  to  be  heavily  associated  to  ion-pairs  and  triple 
ions.  A  theoretical  expression  for  the  triple  Ion  association  constant  similar  to 
the  BJerrum  one  for  ion-pairs  and  to  the  Maaser-BJerrum  theory  of  dimers  has 
been  developed  and  applied  to  the  present  conductance  data.  Ultrasonic  relaxa¬ 
tion  absorption  data  at  much  higher  concentrations  than  the  conductance  data 
reveal  association  to  dimers.*  This  is  also  evident  from  the  microwave  dielectric 
data  showing  no  electrolyte  dielectric  effect  on  the  solvent  and  apparent  lack  of 
presence  of  dipoles.  The  ultrasonic  data  can  be  rationalized  by  a  two-step 
dimerization  mechanism. 
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Introduction 

A  previous  conductance  study  of  the  Ionic  association  of  LIBF4  In  2- 
Methyltetrahydrofuran1  and  in  1,2-Dimethoxye  thane2  of  respective  static  per¬ 
mittivities  e=  6.2  and  €=7.0  has  been  reported.  In  the  same  works, 1,2  the 
diffusions]  rotational  relaxation  dynamics  of  ion  pairs  was  studied  by  dielectric 
relaxation.  In  addition,  ultrasonic  relaxation  revealed  some  dimerization  of  the 
ion-pairs  and  a  kinetic  Investigation  of  this  additional  process  was  reported. 

It  was  of  interest  to  extend  the  above  studies  in  a  medium  of  lower  per¬ 
mittivity  as  D Imethoxymethane  (DMM)  of  e=  2.76  at  25°C,  where  presumably, 
the  electrolyte  exists  completely  as  ion-pairs  at  all  finite  concentrations,  but 
where  the  extent  of  triple  Ion  formation  and  quadruple  formation  Is  sizable  and 
may  become  preponderant  at  high  concentrations  (00.1,  0.2M).  On  the  prac¬ 
tical  side,  LiBF4  dissolved  in  ethers  Is  a  system  relevant  to  secondary  batteries 
construction.  It  was  of  interest  to  report  a  quantitative  study  in  a  solvent  of 
very  low  permittivity  as  DMM  («=2.76  at  25°C)  which  may  constitute  an 
extreme  limit  for  an  electrolyte  solution.  In  the  process,  we  have  derived  a 
new  BJerrum-like  expression  for  the  association  to  triple  ions  as  reported  below. 
This  expression  and  the  corresponding  one  for  dimer  formation,  already 
presented,2  may  give  some  theoretical  guidelines  to  the  extent  of  association  up 
to  the  quadrupoles,  but  below  the  larger  aggregates  which  may  precede  the 
eventual  separation  of  the  electrolyte  from  the  liquid  phase.  For  the  sake  of 
clarity,  after  the  experimental  part  the  conductance  and  theoretical  aspects  lead¬ 
ing  to  the  triple  ion  formation  constant  will  be  dealt  with  first,  followed  by  the 
relaxation  dynamics  study  by  ultrasonic  and  microwave  dielectric  relaxation 
which  involves  itself  with  dimer  formation. 


The  equipment  for  the  conductance15  and  microwave  dieletrlc  relaxation 
has  already  been  described.  For  the  ultrasonic  work  the  pulse  Instrumentation 
has  been  automated,  In  data  capturing,  by  mounting  over  the  dual  crystals 
Interferometric  cell  a  Mltutoyo  Series  164  digital  micrometer  (resolution 
±0.00005  inch)  and  associated  interfaced  digital  counter  and  printer,  giving  a 
hardcopy  of  the  displacements  for  attenuation  increments  (expressed  in  deci¬ 
bels)  of  the  standard  comparison  signal. 

For  the  chemicals,  the  solvent  DMM  (Aldrich)  was  refluxed  over  sodium 
and  benzophenone  until  a  bluish  coloration  Indicating  absence  of  peroxides  was 
present.  It  was  then  distilled  in  the  same  all  pyrex  apparatus  and  used  shortly 
afterwards.  LiBF4  (Aldrich)  was  dried  in  in  vacuo  at  ~60°C  overnight. 


Results  and  Discussion 

a)  Electrical  Conductance 

Figure  1A  reports  the  electrical  conductance  data  In  the  form  of  log10A,  vs. 
logjoc  at  t=25.00°C  for  LiBF4  in  DMM.  Several  runs  with  independently 
prepared  stock  solutions  and  solvent  were  used  In  order  to  insure  reproducibil¬ 
ity  of  the  results.  Figure  IB  reports  the  same  data  for  C<0.1M  elaborated  In 
accord  to  the  Fuoss-Kraus  triple  ions  theory.5 

An  Aq'Kj 


r~  Aq 

A  g(  c)  vc  =  -  ==  + 


v/i<r '  vki  a„)c 


(i) 


where  g(c)  is  a  term  lumping  together  all  the  interionic  terms,  KA  and  I<T  are 


the  ion- pairs  and  triple  ions  formation  constants,  A0  and  Aq  the  limiting  con¬ 


ductivities  of  the  single  ions  and  triple  ions  respectively.  In  the  above,  the  arbi¬ 


trary  condition  Aq  =  — A0  has  been  retained  as  done  previously.1,2  Further, 

3 


for  the  calculations  of  the  Walden  rule  has  been  used  with  the  following 


data  from  the  literature: 


-  In  THF  at  25°C  (viscosity  tj=0.(XM6  poise)6  Xfj=30.6n  ”  'cm^eq'  \  7  hence 
Xf>li7=0.i68. 

In  nitrobenzene  ( 17=0.01823  poise)*  X^=22.m“I  cm*eq"1;*6  hence 
XhV?=0.403. 

It  results  that  In  DMM  ( 17=0.00315  poise)9  X°i=53.3n "‘cm^q"1  and 
X^=127.on"1cm,eq_,t  or,  A£jBp4=l81fl _,cm9eq_1  and  A  J=-|-  A0=12i 
0-1cm2eq_1  In  DMM  at  25°C. 

From  Figure  IB,  the  solid  line  calculated  by  linear  regressions  gives 
r^O.98,  Intercept^ — -  7 .20x10" 7  and  slope  =  =  3.08x10"'*,  from 

v/kT 

which  Ka=6.3x1018M-1  and  KT=828M“1. 

By  equating  KA  to  the  Fuoss-Jagodzlnskl  expression 


Kfj=  — Ld!  e"1/a  eb 
3000 


(H) 


tt  results  In  a  charge  to  charge  separation  In  the  pair  d  =  5.1x10"  ®cm.  Similarly 
by  equating  KT  to  the  Fuoss-Jagodzlnskl10  triple  ion  theoretical  expression 


kfj  _  ,yLa3  e-3/4  ee*/2«kT 
T  1000 


(HI) 


It  results  In  Ion  to  ion-pair  separation  a=  16.9x10"  ®cm,  a  value  that  seems  too 
high  to  be  reconciled  with  the  triple  Ions  model  and  definition  envisaging  the 
three  Ions  at  contact,  whereas,  the  value  a~l.5d=7.7xl0"*cm  would  appear  to 
be  a  more  reasonable  parameter.  Therefore,  we  have  decided  to  re-examine 
the  triple  Ion  theory  for  KT  without  the  constriction  that  the  complexing  ion  be 
In  contact  with  the  pair  to  be  defined  as  a  triple  ion.  We  shall  follow  the  guide¬ 
lines  of  the  BJerrum  theory  for  comparison  and  the  one  for  the  dimers 
(Maaser-BJerrum  theory)  already  presented.9  In  the  past,  dealing  with  the  three 
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ion  interactions  Fuoss  and  Kraus11  have  indeed  derived  a  triple  ion  expression 
based  on  the  BJcrrum  model.  Their  derivation  lead  to  a  integral  which,  to  date, 
is  dim  cult  to  solve,  short  of  resolving  to  a  graphical  integration.  We  thought  it 
worthwhile  to  try  to  arrive  at  a  more  manageable  solution  of  the  problem  by  a 
simplified  derivation  based  on  the  interaction  between  an  ion  and  a  permanent 
dipole  of  moment  p. 

b)  Ion- dipole  interaction 

We  will  start  by  defining  r_  as  the  distance  between  a  given  positive  ion 
and  the  negative  end  of  a  dipole  ion-pair.  Let  r  be  the  distance  between  the 
same  Ion  and  the  center  of  the  dipole  and  r+  the  distance  between  the  Ion  and 
the  positive  end  of  the  dipole  ion  pair.  V  is  the  angle  between  r_  and  the  line 
passing  through  the  dipole  axis,  and  4  Is  the  separation  of  the  two  charges  in 
the  dipole.  We  have  then  (Fig..2A)  r_  -  r+=dcos0\  and  the  ion-dipole  poten¬ 
tial  will  be 

^  _  ef  _1 _ 1  I  _  e_  r+  ~  r-  ~  _  edcosfl' 

c  l  r.  r+ )  e  r+r_  er* 

In  the  above,  e  is  the  static  permittivity  of  the  solvent  and  the  approximation 
r+r_~r2.  If  concomitantly  one  writes  where  0  Is  the  angle  between  r  and 
the  line  passing  through  the  dipole,  one  has 

*3?-££25l  (iv) 

era 

where  p  =  ed  Is  the  dipole  moment  of  the  dipole  separated  by  a  distance  r>d 
from  the  positive  ion.  The  ion-dipole  potential  energy  will  then  be 

(V) 

«ra 

The  probability  of  finding  an  ion  around  a  dipole  taken  as  the  reference  of  a 
coordinate  system,  will  be  a  function  of  the  concentration  of  the  ions 


(Nco/l000)e_u/kT  surrounding  the  dipole;*  of  the  volume  shell  element 
4jrr2dr;  and  of  the  solid  angle  ratio  dw/47rf  expressing  the  point-to-point  spatial 
orientational  probability  between  the  Ion  and  the  central  dipole  (Fig.  2B).  In 
polar  coordinates  r,  0  and  we  will  have 

dw  _  rslnfldtfrrdfl  _  slnfldfld^ 

4*  ”  47rr2  ~  4tt 

and 

(tjieosl  ^ 

«■“  ](«rs<lr) 

The  above  expression  can  be  Integrated  for  t/’—O  to  27T  obtaining  the  probability 
that  the  ion  is  contained  In  a  solid  annulus 


dP(r,0)  =  e«#*eo«#/<r*kT  (4jt  r2dr)  -s-°gdg 

1CXX)  v  '  2 


dP(r.0)  = 


1000 
2ir  Nca 
1000 


i^dr  ( e«pcostf/«r*kT  )  sin^d^ 


Call  cosfl=y,  dy=-slntfd0.  Then,  given 

W 

27rNco 


dP(r)  = 
we  shall  have: 
dP(r)  = 

The  Integral 
t 


1000 


2?rNcg 

1000 


r*drf  (  e*P«»*/«*T  )  sln0d0  t 


o 

i 


r*drj  e(«f‘Ar*kT)ydy  . 


-l 


f  e(fie/cr*kT)ydy  _  1  f  JLut/ir1  kT)y  d(  eA  _ 

J  •  (#ie/crskT)  er2kT 


-i 


Then 


dP(r)  = 


e^eArhcT.e-Me/crhcT  =  asinhlpe/cr^kT] 
(/ie/cr2kT)  (#ie/cr2kT) 


2jt Nc<r  ,  .  [  e"  V«r*kT_  e’ >"'«T,kT 


r*dr 


(VI) 


(VII) 


1000  l  (e/i/erVkT) 

In  order  to  integrate  with  respect  to  r  one  has  to  set  a  finite  upper  limit  to  r. 


•  O  I*  U»«  d«|rt«  of  dtwodaUoa  of  U«  tkclroljrt* 


The  above  function,  for  a  given  finite  thickness  Ar  shows  a  minimum  with  r, 


defined  by 


d  |  AP( r)  |  _  47tNco  1  d  «  (  /ie  |  1  _ 

dr  1  At  j  1000  pe/ekT  dr  rS"  Ui^kT  J  J 


which,  operating  the  derivative  leads  to 


4r3sinh(  ~ — —  +  r * 

i  trakT 


erakT 


2  ne 
er3kT 


4r®  tanhf  — ^ —  |  -  JtL-  r3 = 0 
l  <r2kT  J  cr2kT 

Call  Y=(/ie/cr2kT).  Then  2tanhY-Y=0 


tanhY 


=  0.5 


(VIII) 


which  Is  satisfied  for  Y=1.915«2.  The  probability  function  will  have  a 
minimum  for  Y«2  which  means  physically  that  the  probability  will  go  through 
a  minimum  for  a  squared  distance  r2=q2  that  corresponds  to  an  aligned  (0=0) 
lon-dipole  configuration  of  energy  equal  to  2kT  or; 

-i*r .  <dc) 

For  an  electrolyte  as  LiBF4  of  ji=16xlO“18esu  cm  in  the  solvent  1,2- 
Dimethoxye thane  (c=7.0  at  25°C)  q=11.55xl0~18cm,  whereas  in  DMM 

(e=2.76  at25°C),  by  retaining  the  same  n.  It  results  in  q=18.4xl0_8cm.  Fol¬ 
lowing  BJerrum,  we  shall  integrate  dP(r)  for  r  between  a  and  q  -  a  being  the 
minimum  distance  accessible  to  an  ion  to  approach  a  dipole. 

,  4£NC0  J  2 


1000 


(e/i/e  r*kT) 


and  given 


Y  =  JLH  ...  r  = 
erakT 


\l/a 

_£!_  I  -y—  1/2 
ekT  I 


=  -f  — ] 

{  ekT  J 


Then 


caakT 


dY 


*  /  .3/2 

p/y\  _  _  **Nccr  f  [  1  JL  y-5/8  slnhY 

V  *  ~  1000  Jb  I  ckT  )  2  Y 

b  .  .3/2 

_  2irNc<y  f  |  *ie  ]  slnhY 

P(Y) "  1S55~ •(  I Ht  I 

and  by  multiplying  and  dividing  by  a3 

pm  =  ^Nega3  b3/2  r  slnhY 

1000  J2  y7/a 

The  integral  can  be  solved  by  expanding  in  series  the  convergent  function 

v3  \r5 

slnhY  =  Y +  —  + —  + - 
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(X) 


Q  =  /  -Slnff  dY  =  /  [  Y*s/2  +  —  Y*1/2  +  —  Y3/a  +  —  Y7/*  |  dY 

Ja  Y7^2  j  l  3!  5!  7!  j 

_  Y~3/2  y1/2  Y3^ 

Q  ~  11(3/2)  +  31(1/2)  +  51(5/2)  +  71(0/2)  +” 

b 


Y<a-(S/2)) 

q  -  yi  — i - 

r  (»-(•/*))»* 


.for  all  odd  n's. 


Therefore 


P(Y)  =  2*Nc<ra3  b3/2  Q  =  27rNcqa3  b3/2 


„  Y^n-5^ 

».*a  (n-5/2)n! 


Y*b 


(XI) 


Y-S 


1000  1000 
The  function  Q  Is  a  converging  function  by  increasing  n.  By  limiting  n  to 

n=  15  one  obtains 


b-3/2  bi/a  b«/» 

Q  _  _  |  _  |  _ 

n-ts  1.5  0.5X3!  2.5x5! 


b*7/8 


b21/2 


4.5X7! 
bM/« 


b**  b7 


7x9! 
-  0.2556 


8.5X11!  10.5X13!  12.5X15! 

and  we  shall  identify  the  ratio  P(Y)/ca  as  the  triple  ion  formation  constant  Kx, 


associated  to  either  of  the  two  processes. 
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B8 


AB  +  B"^  ABj 
AB  +  A^  A2B+ 

taken  to  occur  to  the  same  extent  (symmetrical  approximation) 


ar 

^  * 


T  (l-cr-3ax)cac  ca 

namely,  a  degree  of  triple  ion  formation  o T  per  unit  concentration  of  free  ions. 


Then: 


Kt  =  — N—  b3^Q  =  K0b3/aQ  . 
T  1000 


(XII) 


Table  I  reports  the  calculated  values  of  Kt  as  a  function  of  a  for  LiBF4  in  DME 


(=1,2  Dlmethoxyethane)  of  permittivity  e—  7.0  for  which 


H =15.8x10“  18esu  cm.  2  It  was  found2  that  Kt(exp)=50  M~ 1  from  conductance 


experiments.  Judging  from  the  above  figure  it  would  result  su^5.7xi0  'em,  a 


very  reasonable  figure  when  compared  with  the  charge  separation  of  the  dipole 


d=/i/e=(  15.8x10"  18/4.8x10"  l0) =3.3x10"  8cm 


axiom" 


ar=l.5d=4.8xl0“3cm.  The  last  model  envisages  the  triple  ion  as  composed  of 


three  identical  spheres  in  an  aligned  configuration  (corresponding  to  a 


minimum  in  the  potential  energy).  Very  reliable  data  for  KT  as  a  function  of 


the  permittivity  exist  in  the  literature  for  the  system  I-Am4N-N03  in  H20 


dioxane  at  25°C. 


Table  II  compares  the  experimental  data  for  iog10KT  with  the  calculated 


logKT's  and  the  corresponding  values  of  a  used.  The  calculated  KT  corresponds 
to  the  best  fit  by  varying  a.  The  values  of  K^a)  using  the  averaged  a’s  are  also 


shown.  Figure  3  reports  the  log10KT  ys*  1A  for  the  above  system.  The  solid 
lines  are  the  calculated  values  of  log10KT(6>64)>  The  fit  to  the  experimental 


data  is  comparable  with  the  one  obtainable  with  the  Fuoss-Jagodzinski  function 


Kfj  and  the  parameter  a=10.7xl0~*cm 


K»  =  „p(W 

T  1000 


with  /?=e2/2«akT. 

In  fact,  the  average  lA^logK^calc)-  loglC^exp)  |=0.14  with 
a=9.64xl0"8cm,  whereas  |A|=logI<T,-logKT( exp)  {=0.15  with  a=l 0.7x1 0~8cm. 
The  present  approach  however,  does  not  contain  constraints  in  defining  a  triple 
ion  only  when  the  ion  is  "in  contact"  with  a  dipole.  Rather,  it  defines  a  triple 
for  r  varying  between  a  and  q. 

We  then  wish  to  compare  Kx  (Eq.  XII)  with  the  experimental  result  for 
LiBF4  in  DMM  at  25°C.  Using  the  value  of  ^=16xl0*18esu  cm,  obtained  in 
DME  solvent  by  dielectric  relaxation,  Table  I  reports  the  values  of  KT  as  a 
function  of  the  parameter  a.  The  fit  with  the  experimental  value  KT=826M~ 1 
is  obtained  for  a=8.0xl0~8cm.  This  value  is  remarkably  close  to  the  one  calcul¬ 
able  from  the  axiom  a=1.5d=7.7xl<T8cm,  d=5.lxl0~8cm  being  calculable 
from  the  KA  obtained  from  the  conductance  data. 
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for  LiBF. 


Table  I 

Calculated  KT  for  L1BF4  In  DME  (c=  7.0)  and  In  DMM  (c=  2.76)  at 
208.2K  as  a  function  of  the  separation  distance  a  between  ions  and  dipoles. 

Solvent:  DME;  q=(— )I/2=11.47xlO_8cm;  /i=15.8xl0"l8esu  cm 
v  2ckT 


K^exp) 

M~ 1 


axlO8 

b 

K. 

Q 

KT 

1 

(cm) 

— 

(M-1) 

— 

(M"‘) 

1 

a 

4 

16.455 

0.2421 

398.5 

6.4X1 03 

k 

c 

5 

10.531 

0.4728 

12.071 

195. 

5.5 

8.704 

0.6293 

4.115 

66.5 

• 

5.6 

8.395 

0.6643 

3.465 

56.0 

a 

-  5.7 

8.103 

0.7005 

2.955 

47.7 

1 

6 

7.313 

0.8170 

1.463 

31.7 

P 

-♦ 

• 

7 

5.373 

1.2974 

0.828 

13.4 

• 

t 

9 

3.250 

2.7574 

0.3055 

4.94 

m 

<N 

11 

2.176 

5.0345 

0.0532 

.  0.86 

j  ' 

Solvent: 

DMM;  q= 

18.39x10" 

18cm;  £1=16x10" 

18esu  cm 

l 

7 

13.802 

1.2974 

86.34 

5474. 

i 

7.99 

10.594 

1.9294 

12.535 

834. 

8 

10.567 

1.9366 

12.336 

821. 

1 

10 

6.7830 

3.7825 

1.5071 

100. 

12 

4.6965 

6.5362 

0.6283 

41.8 

> 

14 

3.4505 

10.379 

0.3469 

23.1 

16 

2.6418 

15.493 

0.1722 

11.5 

t 

18 

2.0873 

22.060 

0.0268 

1.78 

Table  II 


Experimental  values  of  Iog10Kx  and  calculated  values  of  tog^IC^a)*  for  the  sys¬ 
tem  Isoamylammonium  nitrate  In  H20-dloxane  mixtures  at  25°C.  The  values 
of  logjoK^a)  have  been  approximated  to  logjoK^exp)  by  varying  a  in  steps  of 
0.1x10“  8cm.  The  value  of  k^oKjla)  corresponds  to  the  calculated  Kx  for  the 
averaged  a  for  the  system,  a=9.64xl0_8cm 


€ 

IoSioKt^P) 

lOgloK^a) 

axlO8 

(cm) 

log10KT(a) 

0.60 

2.38 

4.68 

4.65 

9.6 

4.62 

1.24 

2.56 

4.12 

4.11 

9.9 

4.30 

2.35 

2.90 

3.50 

3.50 

10.1 

3.79 

4.01 

3.48 

3.00 

2.98 

9.9 

3.12 

6.37 

4.42 

2.50 

2.51 

9.4 

2.41 

9.50 

5.84 

2.00 

2.01 

8.9 

1.80 

*  The  value  #4=29x10“  J8esu  cm,  corresponding  to  d=6.04xl0“8cm  (#i=de)  has 
been  used  to  calculate  K^a).  a  Is  the  experimental  value  for  the  ion-dipole 
separation  distance,  d  is  the  dipole  charge  to  charge  separation  distance. 
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c)  Dielectric  and  ultrasonic  relaxation 

Figure  4A  reports  the  real  part  t'  of  the  complex  permittivity  e*  =  t'~  Je" 
plotted  vs.  the  frequency  f  for  L1BF4  0.35M  In  DMM  at  25°C.  The  solid  line 
Is  the  fitted  function  which  appears  capable  to  interpret  the  available  data 
according  to  a  single  Debye  relaxation  process 

t'  =  (£»-‘“)7^  +  e~  •  (xra) 

Figure  4B  is  the  Cole-Cole  plot  of  e"  vs.  e1  and  the  semi-circle  corresponds  to  a 
single  relaxation  process  with  e'  given  by  the  equation  above  and  ew 


(«o-  £oo) 


(XIV) 


l+(fAr)2 

with  parameters  eQ  =  2.78,  =  2.02  and  tT  =  60  GHz.  Notice  that  Saar  et.al.® 

reported  for  the  solvent  DMM  c0  =  2.76  =  2.18  and  fr  =  75  GHz.  Hence, 

the  presence  of  0.35M  Is  only  detectable  in  a  shift  of  the  solvent  relaxation  to 
lower  frequencies  and  a  change  of  the  relaxation  strength  (e0-  c^,)  which  is 
apparently  due  mostly  to  a  change  in  The  latter  one  is  only  an  extrapolated 
parameter  because  of  possible  effects  op  non-Debye  nature  at  f>100  GHz.12 
The  rather  remarkable  findings  for  the  present  system  is  the  absence  of  a 
"solute"  relaxation  around  1-3  GHz  which  has  been  the  observed  behavior  for 
alkali  salts  in  ethereal  solvents  studied  so  far.13,1’2  This  apparent  Invisibility  of 
the  dipolar  pairs  L1BF4,  despite  the  very  large  KA  found  by  conductance,  can  be 
interpreted  in  two  ways:  either  the  pairs  are  so  heavily  solvated  that  their 
diffuslonal  rotational  relaxation  frequency  is  below  the  range  accessible  to  our 
microwave  measurements,  or  LiBF4  is  heavily  dimerized  above  — 0.1M  to  apo- 
lar  or  antiparallel  dimers  which  are  practically  invisible  by  dieletrlc  relaxation. 


Consideration  of  the  second  hypothesis,  namely  a  sizeable  dimerization,  is 
suggested  by  the  fact  that  already  In  DME2  (e  =  7.0)  an  ultrasonic  relaxation 
was  interpreted  as  due  to  dimerization.  In  DMM  solvent  (e  =  2.76)®  simple 


electrostatic  considerations  would  suggest  the  quadrupole  or  dimer  formation 
constant  to  be  much  larger  than  In  DME. 


Figure  5  shows  a  representative  plot  of  the  excess  sound  absorption  per 
wavelength  fi  =  <*excX  plotted  vs.  f  for  L1BF4  in  DMM  at  25°C.  Although  the 
data  extends  to  about  3  decades  In  frequency,  a  single  Debye  process  (solid 
line)  appears  adequate  to  represent  the  data  according  to  the  function 


M  max 


f/fr 


(XV) 


l+CW2 

where  //  =  nmxt  at  the  relaxation  frequency  f  =  fr,  /i  =  aexcX  =  (a-Bf2)u/f, 
and  a  Is  the  sound  attenuation  coefficient,  X  the  wavelength  X  =  u/f,  n  is  the 
sound  velocity  and  B  =  ( a/f2)f>>(r  the  background  (a/f2)  ratio  for  frequencies 
f»fr.  Table  III  reports  the  relaxation  parameters  Hmax,  fr,  B  and  the  sound 
velocity  u  for  all  the  solutions  and  temperatures  Investigated. 


From  Table  III  it  is  apparent  that,  within  experimental  error,  fr  Is  indepen¬ 
dent  of  concentration  at  25°C,  but  that  /im  is  not  linear  with  concentration. 
This  last  observation  negates  the  interpretation  by  a  scheme  A  B  which 
could  represent  a  first  order  (or  pseudo  first  order)  lntra-molecular  process 
(involving  the  solvent). 


Rather,  we  propose  a  multistep  dimerization  process  of  the  type 


kx  k2 

2M  ^  M  •  •  •  M  M2  (XVI) 

k_,  k_2 

with  M  =  LiBF4,  the  monomer  Ion-pair.  This  scheme  leads14  to  two  relaxation 
times  that  for  kx,k_  x»k2,k_2  read 

rf 1  =  4kx(M)  -I-  k_x 

Ti'’  =  + ■ with  K- = (k-/k-)°ii^rl(xvn) 


For  4(M)»K_X 


Tii  =  k-2  +  k2  =  k_2(l+K2)  ,  or 


t,7  1  =  e+“-’‘«'  e-4H-’S/RT(i+Ka) 

h 


(XVIII) 


with  I<2  —  k2A-2*  The  above  Implies  that  a  plot  of  lntr^/T)  vs.  1/T  will  have 


a  slope15 

dlnOV/T) 

d(l/T) 

withK2  =  e-AHa/BTeAS*/R. 


AH  f2  K2  AHs 
R  l+K,  R 


.  (XIX) 


Figure  8A  shows  a  plot  of  ln(rj^  */T)  vs.  1/T.  The  solid  line  was  calculated 
by  linear  regression  giving  r2  =  0.95,  Int  =  16.01  and  Slope  =  -468  from 
which  ASj^  =  -15.4  cal/K  mol  and 


-  486  = 


AH4  Ka  AHS 
R  1+K2  IT 


(XX) 


Also,  Eq.  (XVIII)  gives 


5.68.10'  =  e4sS/B  (l+K,)  .  or 


Table  III 

Ultrasonic  parameters  pm,  fr,  B  and  sound  velocity  u  for  all  the  concentra¬ 
tions  of  LiBF4  in  DMM  and  temperatures  investigated. 


t 

C 

/*roc 

fr 

BxlO17 

UX105 

(oC) 

(M) 

X10S 

(MHz) 

(cm_1s2) 

(cms*1) 

25 

0.35 

430 

85 

48 

1.069 

25 

0.24 

450 

80 

42 

1.065 

25 

0.108 

380 

80 

38 

1.064 

25 

0.093 

280 

80 

36 

1.063 

15 

0.25 

410 

80 

40 

1.142 

5 

0.25 

470 

75 

36 

1.153 

-1 

0.25 

480 

70  . 

38 

1.179 

*uxi05= 1.179-0.0041,  (r2  =  0.8),  for  L1BF4  0.25M  between  25°C  and  -1°C 
reproduces  the  experimental  data  with  a  %  error  =  ■  ‘  U-i  %  =  ±  1 


For  a  two-step  dimerization  process  (XVI),  the  maximum  sound  absorption  per 
wavelength  for  the  slow  relaxation  process,  is  bound  to  the  concentration  of 


the  various  species  by  the  relation14 


-  n  r. 

“  20s  RT  11 


(XXII) 


with  0S  =  (/> u2)-  1  the  adiabatic  compressibility,  p  the  density,  AVS1,  the  isoen- 


tropic  volume  change 


AVSI,  =  AV2  +  - - AV. 

SI1  2  1/4(M)  +  (M...M)  1 

with  AV j  and  AV2  the  isoentropic  volume  changes  associated  with  the  steps  of 


process  (XVI).  Also,  the  function  rn  is:14 


r"  -  [  (Mj  +  W)  +  (M-'-M)  j  (XXII,) 

(the  factor  1/2  Instead  of  1/4  in  front  of  (M)  in  Eqs.  (XXII)  and  (XXIII) 

reported  in  reference  14  was  a  misprint).  If  (M)<(M...M),  (M2),  as  the 

dielectric  data  seem  to  imply,  Eq.  (XXIII)  can  be  approximated  to 


r„  = 


with  K2  = 


(Mg) 

(M...M) 


— i - h  -  -1  -  T‘  = 

(M2)  (M...M) 


l+Ka 


(M...M) 


■.  On  the  other  hand,  one  can  write 


(M...M)  = 


(M...M)  +  (M2) 
1+K2 


(XXIV) 


K2 

r„  s - 2 —  c11"”"  (xxiv) 

(i+k,)! 

where  cd,mer  =  (M...M)  +  (M2)  (Notice  also  that  AVSI,  is  concentration  depen¬ 


dent  and  that  indeed  fin  is  not  linear  with  c). 
Introducing  XXIV  into  XXII 


it  (AVS„)!  K,  „ 

fill  =  — - - c 

‘  20s  RT  (l-fIC2)2 


'.;v 


37  - 


ln(/i„T/u2)  =  In 


IT  pc 


dimer^y2 


Sll 


2R 


+  InK^-  2In(  1 +I<2) 


By  neglecting  the  temperature  dependence  of  the  quantity  pcdlmerAV|n,  the 
above  implies  that 

dln(ji„T/u2)  AH2  2Ka  AH-  AH-  IC-- 1 

- .L  g - 1  + - 1 - L  - - -  — -  (XXVI) 

d(l/T)  R  1+Ka  R  R  Kj+1 

namely,  the  slope  of  a  plot  of  ln(puT/u2)  vs.  l/T  should  be  expressed  by  Eq. 
(XXVI).  Figure  6B  shows  this  plot.  Linear  regression  gives  Intercept  = 
-21.05,  Slope  =  >526.0.  Therefore, 

AH-  K--1 

-526.0  =  - -  — -  (XXVII) 

R  Kj+1 

Eqs.  (XX),  (XXI)  and  (XXVII)  can  now  be  correlated.  Trial  and  error  choice15 
of  1<2  leads  to  compatibility  for  =  7.8  which  leads  to  AH2  =  -1.45K  cal/mol 
AH  =  2.2K  cal/mol. 

Table  IV  reports  the  above  and  the  derived  thermodynamic  and  kinetic 
parameters  using  the  relations  AG2  *=  -  RT  In  =  -  1.22K  cal  at  T  = 
208.2K,  AS2  =  (l/T)(AH2-  AG2)  =  AS^-  ASJ?£  and  AH2  =  AH^-  AH  f2. 

It  would  appear  desirable  to  evaluate  Kq,  the  overall  quadrupole  or  dimer 
formation  constant.  From  scheme  (XVI)  and  the  definitions 

(M2) 


_  (M...M)  _ 

1  (M)2  *  ^  “  (M...M) 

(M...M)  +  (M2) 


and 


ke  = 


it  results  that 


(M)2 

Ke  =  K.d+Ka)  '  (XXVIII) 

and  we  have  obtained  I^  =  7.8.  K,  could  be  evaluated  in  principle  from  the 
Maaser-BJerrum2  theoretical  expression 


—  KoqbqQq 


(XXIX) 


with 


^oq 


47rLa8 

3000 


ea*kT 


and 


.  38  - 


Q,  ■  0.8067  4  +  (b"-(l.5)n)  ,  (XXX) 

^  bq  |  (n+2)!n 

for  odd  n*s. 

For  c  =  2.76  and  by  retaining  ft  =  16x10“  18esu  cm,  we  have  calculated  Q 
extending  n  to  n  =  17.  The  calculation  of  KqBj  is  shown  on  Table  V.  Unfor¬ 
tunately,  KqBj  appears  to  be  too  sensitive  to  the  choice. 


Table  IV 

Thermodynamic  and  kinetic  parameters  for  the  slow  relaxation  dimerization 
process  for  LiBF4  In  DMM 


7.8 

-1.22  K  cal /mol 
-1.45  K  cal/mol 
-0.8  cal/  K,  mol 
2.2  K  cal/mol 
-15.4  cal/K,  mol 
0.76  K  cal/K  mol 
-16.2  cal/K,  mol 
5xl08s“ 1 


of  the  parameter  Sq,  the  minimum  dipole-dipole  separation  distance  to  make  an 
assignment  of  KqBj  s=.K,  short  of  being  an  arbitrary  assumption. 


Table  V 

Calculated  values  of  Kq^  according  to  the  expression  K^  —  K^bQq  as  a 
function  of  aq  the  dipole-dipole  minimum  separation  distance  for 
ft  —  16x10“  ,8cm  esu  the  dipole  moment  of  the  pairs,  T  *=  208.2,  e  =22.76  and 

bq  ts  n*/ea*kT.  ICqBj  =  0  for  a  separation  distance  q  *=  (—■  “~r)1/3  = 

3  CK1 

11.45xl0“  cm,  for  the  above  parameters. 


V 

Qq 

KqBJ 

5 

18.035 

0.3152 

5507 

3.1x10* 

6 

10.437 

0.5447 

23.34 

132.7 

7 

6.572 

0.8640 

2.013 

16.57 

8 

4.403 

1.201 

1.231 

7.00 

0 

3.002 

1.830 

0.6035 

3.04 

10 

2.254 

2.522 

0.3734 

2.12 

11 

1.604 

3.356 

0.1120 

0.642 

Conclusion 

The  Fuoss-Kraus  triple  Ion  conductance  theory  allows  for  the  determina¬ 
tion  of  I<A  =  fixlO^M"1  and  KT  =  826M-1,  pending  the  validity  of  the  Walden 

-  o 

rule  and  the  assumption  A.J  —  — Ae.  The  Fuoss-Jagodzinski  association  rela- 

3 

tlon  for  I<A  gives  d  =  5.1x10"  8cm  for  the  charge  separation  of  the  Ion- pair. 
The  theoretical  expression  XII  derived  above  for  ICT  gives  a  =  8.0xl0“*cm  for 
the  minimum  separation  ion  distance  between  the  ion  and  the  dipole,  close  to 
the  axiom  a  =  l.5d  =  7.7x10“ 8cm.  The  dielectric  data  shows  absence  of  a 
dielectric  relaxation  for  the  solute  in  all  the  frequency  ranges  Investigated.  This 
could  be  rationalized  by  either  the  dipole  pairs  relaxing  below  l  GHz,  or,  more 
likely,  by  the  electrolyte  being  dimerized  to  apolar  dimers. 

An  ultrasonic  relaxation  of  Debye  type  is  interpreted  as  due  to  the  second 
step  of  a  coupled  two-step  dimerization  process.  Kinetic  and  thermodynamic 
parameters  have  been  extracted  from  the  data.  In  particular,  the  formation 
constant  for  the  second  step  of  dimerization  K 9  «=  7.8.  The  constant  for  the 
first  step  I<j  cannot  be  determined  without  a  reliable  knowledge  of  a^  the  pair 
to  pair  minimum  separation  distance. 
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ADDENDUM 


After  the  completion  of  the  present  work  we  have  performed  static  permit¬ 
tivity  measurements  at  f  =  3.5  MHz  with  a  Bontoon  resonator  and  a  cell  of 
capacity  Ce  =  5.07±  0.07p  Farad  of  a  solution  0.093M  L1BF4  in  DMM  at  25 °C. 
The  average  of  two  experiments  gave  e0  —  2.82d;0.0|,  a  value  close  to  the 
figure  extrapolated  from  the  microwave  range  c0  =  2.78  as  reported  above. 
This  confirms  the  hypothesis  advanced  above  that  the  electrolyte  is  heavily 
associated  to  apolar  dimers. 
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Real  part  of  the  complex 
permittivity  €* *  £'-J €’ 
vs.  frequency  f.  for  LiBF4 
0.35M  in  DMM,  t-25*C 

—  Single  Debye  with  2-7®.  €m-  2.02,  fr  -  60GHz 


Cole-Cole  plot  for  LiBF4 
0.35 M  ktDMM.  t-25*C 

- Single  Debye  process  with 

Co"  2.78.  em-2.02,  fr-60GHz 
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Lamb  plot  for  LiBF4 
0.25M  in  DMM 


